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Abstract 

The general strategy, as well as channel-specific details, applied to the measurement 
of the top quark mass at the Tevatron in Run I are reviewed, and the combination of 
the results obtained by the CDF and D0 collaborations presented. The accelerator 
and detector upgrades for Run II are described, and expected improvements in the 
systematic uncertainties on the measurement evaluated. 



1 Introduction 



The mass of the top quark is a fundamental parameter of the standard model. 
Its precise measurement will not only improve understanding of the model, but in 
combination with other experimental data it will also allow to further check the 
model's consistency. In addition, combined with the W mass measurement, it can 
be used to constrain the mass of the Higgs boson. 

This paper is organized as follows: Section || addresses top quark produc- 
tion and decay topologies at the Tevatron; top-quark mass measurements performed 
in Run I are reviewed in Section |3], and the combined CDF and D0 result is given in 
Section |j. This is followed in Section || by a description of accelerator and detector 
upgrades for Run II, and in Section || by the prospects for reduction of systematic 
uncertainties on the mass measurement. Conclusions are drawn in Section [7]. 

2 Top Quark Production and Decay at the Tevatron 

At the Tevatron, top quarks are produced mainly in pairs. At y/s = 1.8 TeV, about 
90% of the cross section is due to qq annihilation and 10% to gluon fusion. The it 
production cross section has been measured by both the CDF HP and D0 ED collab- 
orations to be a{m top = 175 GeV/c 2 ) = 6.5t\'A pb, and a{m top = 172 GeV/c 2 ) = 
5.9 ± 1.7 pb, respectively. These results are in good agreement with the theoreti- 
cal predictions ®, which range from 4.7 to 5.5 pb. Single top quarks can also be 
produced in electroweak processes, but due the smaller cross sections and larger 
backgrounds, these have not yet been observed. 

Because of its large width, the top quark decays to a b quark and a W 
boson before it hadronizes. The various possible W boson decays determine the 
separate channels for mass measurement: events in which both W bosons decay 
hadronically comprise the "all hadronic" channel, and account for « 44 % of all 
produced it events. Events in which one or both of the W bosons decay to ev or [w 
form the "lepton+jets" (« 30 %) and "dilepton" (« 4.9 %) channels, respectively. 
Events in which one or both W bosons decay to rv are not used in the mass and 
cross section measurements. 

3 Top-Quark Mass Measurements 
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3.1 General Strategy 



While different experiments and different channels require specific analyses, these 
measurements have enough in common so that one can speak of a general strategy 
to extract the mass of the top quark from data. This involves four steps: 

• Apply a set of criteria to select the candidate sample, while being careful to 
choose cuts that introduce little bias in the extracted mass. 

• Fit each event's kinematics to the top quark pair production hypothesis, ex- 
tracting a fitted top-quark mass m/j t . 

• Signal distributions are generated for many values of m top , and a likelihood 
(£) fit is used to compare the measured mf it distribution to summed signal 
and background distributions for each value of m top . The fit determines the 
relative contributions of signal and background. 

• Fit the obtained InC versus m top , and extract the mass at the minimum of 
— InC 

In all cases, ti signal distributions have been generated using the HER- 
WIG @) Monte Carlo generator. The background is dominated by multiple-jet pro- 
duction and the production of a W boson in association with additional jets. The 
former is derived from data while the latter is modeled using the VECBOS D event 
generator. 

3.2 All-Hadronic Channel 

Only CDF has published a mass measurement ® in the all-hadronic channel. Events 
are required to have six or more jets, of which at least one is tagged as a 6-quark 
jet by the detection of a secondary vertex ("SVX tag"). Imposing that the masses 
of the top and anti-top quark be equal, and that the invariant mass of pairs of jets 
forming W bosons equals the P^-boson mass, yields a 3-C fit. However, for each 
event, multiple combinations are possible in the association of jets to W bosons or b 
quarks, so the combination resulting in the lowest \ 2 ls chosen as the right kinematic 
fit. The sample consists of 136 events (corresponding to an integrated luminosity 
of 109 pb^ 1 ), with the background contamination estimated at 108±9 events. The 
measured top mass is m(top) = 186.0 ± 10 (stat) ± 5.7 (syst) GeV/c 2 . 
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3.3 Lepton + Jets Channel 

The D0 lepton + jets analysis 0) requires an isolated central lepton, four or more 
jets and at least 20 GeV missing transverse energy. Since the neutrino's longitudinal 
momentum is not measured, the kinematic fit to the it hypothesis is now a 2C fit. 
D0 also computes a discriminant D, based on 4 variables weakly correlated with the 
top-quark mass. This discriminant is computed through a comparison of signal and 
background using classical methods, as well as a neural network. The likelihood fit 
to determine the top mass is then performed on a two-dimensional grid of D versus 
the top mass. Both discriminants yield results that are in good agreement with each 
other, which when combined give m(top) = 173.3 ± 5.6 (stat) ± 5.5 (syst) GeV/c 2 . 
The background contamination is estimated to be approximatively 50 out of the 77 
events that pass all selection criteria. 

CDF 0) uses a similar procedure, but finds that optimal results can be 
obtained by dividing the sample into 4 subsamples: events with two SVX tags 
(5 events), events with a single SVX tag (15 events), events without SVX tags 
but where one jet is identified as a 6-quark jet through a soft lepton tag ("SLT 
tag") (14 events), and finally events without any tags (42 events). The samples are 
fitted separately and the results, shown in Fig. [IJ, when combined yield m(top) = 
175.9 ± 5.1 (stat) ± 5.3 (syst) GeV/c 2 . 

3.4 Dilepton Channel 

As indicated in Section |2|, the dilepton channel has fairly low statistics. On the other 
hand, the background is small, dominated by Z boson decays to charged leptons and 
particle misidentification. The difficulty arises from the presence of two neutrinos, 
for which only the combined transverse momentum can be measured. The system is 
thus kinematically underconstrained, and both experiments use the so-called "neu- 
trino weighting" technique to measure the top quark mass: for every chosen top 
quark mass, there is a sampling of allowed neutrino rapidities. Given the neutrino 
rapidities, the assumed top mass, and the charged-lepton and b-quark momenta, 
the system can now be solved for the transverse and longitudinal momentum com- 
ponents of the neutrinos (with a two- fold ambiguity for each neutrino). A weight is 
then assigned to each solution based on the agreement between the calculated and 
measured missing transverse momentum. All weights are summed for all possible 
combinations of neutrino rapidities. This is done as a function of top mass value for 
each event. 

D0 D has observed 6 events in a sample corresponding to an integrated 
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Figure 1: CDF lepton + jets analysis: reconstructed mass distributions in each of 
the four subsamples (see text). Each plot shows the data (points), the result of the 
combined fit to top + background (dark shading), and the background component of 
the fit (light shading). The insets show the variation of the log-likelihoods with m top 
for the separate fits. 



luminosity of about 125 pb~ l . For each event, the distribution of weights is split 
into five mass regions, and used to form a vector for calculating the mean mass 
value for all events in the sample. Fig. [| summarizes the results based on these 
events. D0 also evaluates the mass in this sample using an extended version of the 
Dalitz-Goldstein and Kondo method. Both results are compatible, and the combined 
measured top-quark mass is m(top) = 168.4 ± 12.3 (stat) ± 3.6 (syst) GeV/c 2 . 

Eight events are observed by CDF 0) (in a sample corresponding to an 
integrated luminosity of 109 pb^ 1 ). Each event's top mass estimate is determined 
from the peak of the weight distribution, and the resulting distribution compared 
with Monte Carlo simulation, leading to a measured top-quark mass of m(top) = 
167 A ± 10.3 (stat) ± 4.8 (syst) GeV/c 2 . 

3.5 Systematic Uncertainties 

The main contributions to the systematic uncertainties on the Run I top-quark 
mass measurements are given in Table |l| The dominant contribution is due to the 
determination of the jet energy scale. Other important contributions come from 
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Figure 2: D0 dilepton analysis: summed event weight function J2i Wi for the 
data sample (points), the sum of the fitted signal and background (solid), and the 
background alone (dashed). The error bars indicate the rms observed for five-event 
samples in ensemble tests. The inset shows the corresponding fit to —InC, drawn as 
a solid curve in the region considered in the fit. 



uncertainties in initial and final state radiation, and modeling of the background. 
All of these will be discussed in more detail in Section §. 

4 Combined CDF and D0 Measurement 

The Top Averaging Group HID has combined the 5 mass measurements described 
in Section [|. As expected, after taking into account all correlations, the best con- 
tribution comes from the lepton + jets channel. The combined result is m(top) = 
174.3 ±3.2 (stat) ±4.0 (syst) GeV/c 2 , which is compared with the individual mea- 
surements in Fig. H 

5 Run II Upgrades 

Run II is formally scheduled to start at the Tevatron on March 1st, 2001. This 
will follow a shutdown during which both the accelerator complex and the collider 
detectors will have undergone significant upgrade from their Run I configurations. 
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Table 1: Summary of main contributions to systematic uncertainties (in GeV/c 2 ) 
in the measurement of the mass of the top quark in Run I: a) jet energy scale, b) 
model for signal (mainly initial and final- state radiation), c) Monte Carlo generator, 
d) Uranium noise and multiple interactions, e) model for background, and f) method 
for mass fitting. Contributions to the uncertainties differ somewhat between the two 
experiments. For example, for CDF measurements, the uncertainty due to multiple 
interactions is included in the uncertainty in the jet energy scale. 
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5.1 Upgrade of the Accelerator Complex 

In Run II, the Tevatron is expected to operate at a center of mass energy y/s = 
2 TeV, an increase of 200 GeV over Run I. For top quark studies, this yields a 
substantial rise of 40 % in the it cross-section. 

In addition to that, the completion of the Main Injector should provide 
significant improvement in the luminosity. Both experiments are expected to record 
data corresponding to an integrated luminosity of 2 fb' 1 in the first two years, 
followed, after a brief shutdown, by 5 fb^ 1 per year until the start up of the LHC. 

5.2 D0 Upgrade 

The main features of the D0 upgrade involve the addition of a 2 Tesla solenoid 
inside the calorimeter, combined with improved tracking and particle identification 
capabilities. 

Two central tracking systems have been built: an extensive silicon detector 
consisting of barrels for central tracking and disks covering the forward regions, and 
a scintillating fiber tracker between the silicon and the solenoid. The silicon barrels 
have four layers mainly of double-sided silicon sensors, while all disks provide two 
hit projections. The fiber tracker consists of 8 coaxial cylinders, and is used in the 
Level- 1 trigger. 

In addition to the tracking systems, preshower detectors will be used to 
identify electrons and photons. Also, the muon coverage has been extended and a 
new three-level trigger system is being built. 
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Tevatron Top Quark Mass Measurements 
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Figure 3: Individual and combined top-quark mass measurements at the Tevatron. 
5.3 CDF Upgrade 

The upgraded CDF detector will have improved central tracking detectors, as well 
as increased coverage in pseudorapidity for lepton identification. 

The new silicon detector consists of three parts: "Layer-00", silicon bar- 
rels, and the "Intermediate Silicon Layers (ISL)". Layer-00 is a single layer of silicon 
located very close to the beam, leading to significant improvement in impact pa- 
rameter resolution. The barrels have five layers of double-sided silicon sensors, with 
the ISL located at larger radii, providing additional points between the silicon and 
the Central Outer Tracker (COT). The COT is an improved version of CDF's Run 
I central drift chamber: it has 12 layers of wires, less material, and operates with 
shorter drift times. 

New end plug calorimeters have been built, giving better coverage for elec- 
tron identification, and the muon coverage has also been increased. In addition to 
this, a time-of-flight system is located between the COT and the solenoid, and a 
new three-level trigger system is being implemented. 
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6 Systematic Uncertainties on the Top-Quark Mass Expected in Run II 

In this section, we discuss ways to reduce the dominant systematic uncertainties on 
the mesurement of the mass of the top quark. 



6.1 Determination of Jet Energy Scale 

For Run I data, CDF and D0 use a somewhat different approach to determine the 
jet energy scale. In D0, the essential steps of the procedure are the following: the 
electromagnetic scale is normalized by the mass of the Z-boson, and the hadronic 
response is then derived by balancing 7 + jets events. Corrections for the underlying 
event, multiple interactions and noise, are determined from minimum-bias data, 
while corrections for leakage outside the jet cone are computed based on test-beam 
data and Monte Carlo simulation. Finally, the conversion from jet to parton energy 
is based on Monte Carlo (HERWIG), with separate treatment of light and 6-quark 
jets. A small remnant pseudo-rapidity dependence is corrected using data from 7 
+ 1 jet events, and the results are cross-checked using [Z — > ee) + jets events. In 
D0, the required corrections are determined separately and symmetrically for data 
and Monte Carlo, and then the relative difference between data and Monte Carlo 
is added to the systematic uncertainty. Since the top-quark mass is measured by 
comparison of data with simulation, this is a valid approach. 

While many of the steps are similar in CDF, there is one fundamental 
difference: in CDF, the Monte Carlo simulation, and in particular the fragmentation, 
is tuned to the data, exploiting the available accurate tracking information. As is 
shown in Table [TJ, the two different approaches do not seem to lead to significantly 
different results. 

In Run II, the additional and improved data will yield at least two new 
handles on the jet energy scale in it events. The first is given by the larger sample 
of lepton + jets events with two 6-quark tags. In 2 fb" 1 , this sample is expected 
to contain of the order of 300 to 500 events, allowing accurate calibration of light 
quark jets through the reconstruction of the W-boson mass. Probably the most 
precise way to determine the calibration of the jet energy scale is given by the large 
sample of Z — > bb events expected in Run II. This type of events has already been 
observed by CDF in Run 1 03) , but the sample is unfortunately to small to be used 
for determining the jet energy scale. It is expected that using these events in Run 
II will allow reduction in this dominant uncertainty to below 1%. 

6.2 Initial and Final-State Radiation 

A second sizeable contribution to systematic uncertainty is due to the treatment of 

initial and final-state radiation. Clearly, additional jets from initial-state radiation 

should not be included in the kinematic fits, but those due to final-state radiation 

certainly should be. Unfortunately, such additional jets are the least energetic only 
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about 55 % of the time. However, theoretical work is progressing, and present and 
future versions of the Monte Carlo generators should reproduce these processes with 
better accuracy. Furthermore, the sample of lepton + jets events with extra jets 
in Run II should be large enough to allow more detailed study and comparison to 
Monte Carlo. This can be done using variables such as the number of additional 
jets, relative angles between jets, etc. 



6.3 Backgrounds 

Multiple-jet QCD background is determined directly from data, and does not repre- 
sent a major problem. In Run I, a significant systematic uncertainty came from the 
H^+jets background, which was modeled using the VECBOS event generator. In 
the meantime, W and Z-boson transverse momentum distributions at the Tevatron 
have been studied in detail , contributing to our understanding of these vector 
boson production processes. Better Monte Carlo generators already exist, so this 
uncertainty should be reduced appreciably in Run II. 



7 Conclusions 



In Run I at the Tevatron, the top quark mass has been measured with an accuracy 
of 3 %. The measured value is m(top) = 174.3 ± 3.2 (stat) ± 4.0 (syst) GeV/c 2 , 
which corresponds to a top Yukawa coupling Xt op (Q 2 = ^f op ) = 1.00 ±0.03. 

In Run II, thanks to the improved accelerators and detectors, statistics 
will not be an issue. Furthermore, the additional data will provide better control 
of systematic uncertainties, where the main contribution will probably shift to the 
treatment of initial and final-state radiation. A challenging goal will be to reduce 
the systematic uncertainty on the top quark mass to less than 2 GeV/c 2 . 
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